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Abstract 

Vacuum characteristics quantifying dynamical tendency toward self-duality in gauge 
theories could be used to judge the relevance of classical solutions or the viability 
of classically motivated vacuum models. Here we decompose the field strength of 
equilibrium gauge configurations into self-dual and anti-self-dual parts, and apply 
absolute X-distribution method to the resulting polarization dynamics in order to 
construct such characteristics. Using lattice regularization and focusing on pure-glue 
SU(3) gauge theory at zero temperature, we find evidence for positive but very small 
dynamical tendency for self-duality of vacuum in the continuum limit. 

1. Introduction. The availability of topologically non-trivial classical solutions to Yang- 
Mills equations pj, combined with technology of semiclassical calculations, provided the 
framework for discussion of important non-perturbative effects such as the Ua{^) problem 
of QCD [2], 6'-dependence [3] or the baryon number violation in the Standard Model [2]. 
However, the merits of semiclassical approximations have been questioned since the early 
days, especially at zero temperature and away from the Higgs phase [1] . 

The applicability of semiclassical approach is related to vacuum properties of the theory 
in question. Indeed, if classical solutions represent a sensible starting point of a systematic 
analysis, then it is expected that their properties will manifest themselves in configurations 
dominating the associated path integral. In that vein, the semiclassically motivated picture 
of QCD vacuum based on instantons has been evolving [5l E] , focusing mainly on the physics 
of spontaneous chiral symmetry breaking [7] and with confinement being rather problematic. 
In fact, the large N-motivated arguments of Witten |1] qualitatively amount to the propo- 
sition that, in the confining vacuum, the condition of field being pure gauge on a boundary 
enclosing 4-dimensional subvolume cannot be satisfied even approximately. This implies 
that the topological charge of true QCD vacuum is not expected to be locally quantized 
in near-integer values as predicted by instanton-based picture. The absence of this local 
quantization is consistent with direct lattice QCD calculations [8J. Moreover, it has been 
found [9J and confirmed [10] that the topological charge in pure-glue QCD vacuum organizes 
into a lower-dimensional sign-coherent global structure. This signals vacuum arrangement 
of very different nature than the one suggested by semiclassical reasoning. 

It would be convenient and desirable to have a vacuum-based measure at our disposal, 
capable of probing the "degree of classicality" in the theory. What we have in mind is a 
simple characteristic, computable by lattice methods, that could signal the basic expectation 
for validity of semiclassical approximations and the semiclassical picture of the vacuum. 



Possible strategy for such an approach has been put forward in Ref. [TT], but it has not been 
properly carried out. In particular, for gauge theories, classicality is intimately connected to 
(anti) self-duality. Indeed, self-duality of a field configuration suffices to solve the equations 
of motion in pure-glue SU(N) gauge theory (see e.g. \12\ for recent review), and the true 
local minima of the action are of this type. Thus, in practical terms, self-duality of vacuum 
configurations is almost synonymous with what we wish to characterize. 

In their initial study of self-duality, the authors of Ref. [11] pursued an indirect approach 
via Dirac eigenmodes. In particular, they proposed a differential characteristic of local chiral 
behavior in the near-zero modes (X-distribution), which could also be viewed as a low- 
energy reflection of possible self-duality in the gauge fleld. In Ref. [I3] the same technique 
has then been used in conjunction with a particular form of truncated eigenmode expression 
for fleld strength itself. However, the crucial shortcoming of the above attempts is that they 
are arbitrary with respect to the choice of the polarization measure. Indeed, as demonstrated 
in Ref. [11], it is possible to obtain any desired behavior of X-distribution by the appropriate 
adjustment of the reference frame (polarization function). Thus, the method in this form is 
inherently kinematicalQ 

This issue has been recently addressed in detail, and resulted in the general polarization 
method of absolute X-distribution [Tl]. The polarization measure is uniquely flxed in this 
framework, so as to provide a differential comparison of polarization tendencies in ques- 
tion relative to those involving statistically independent polarization components. By virtue 
of being a valid polarization measure and, at the same time, a comparator to statistical 
independence, absolute X-distribution Pa{X) represents a dynamical polarization charac- 
teristic. The associated integrated quantity, the correlation coefficient of polarization Ca 
(— 1 < Ca < 1), expresses the overall dynamical tendency for polarization. 

One noteworthy virtue of adopting correlational characteristics is that it alleviates the 
usual worries about distortions of vacuum structure due to "ultraviolet noise" . These worries 
typically result in attempts to suppress high frequency fluctuations, and were in fact part 
of the reason for focusing on low-lying Dirac modes in the above context of self-duality. 
However, with absolute polarization methods, this is not necessary and it is appropriate to 
approach the problem directly in equilibrium conflgurations of the gauge fleld thus capturing 
the properties of full dynamics. Indeed, if ultraviolet fluctuations represent pure noise, they 
will not affect the direction of dynamical tendencies. On the other hand, if these fluctuations 
exhibit correlations, then there is no reason to suppress them since they can affect the physics, 
even at low energy. 

Given the above, we propose the correlation coefficient of polarization Ca (and the un- 
derlying absolute X-distribution) , evaluated for the duality decomposition of the gauge fleld 
strength tensor, as a suitable vacuum measure of self-duality, and thus a candidate for com- 
putationally accessible characteristic of classicality in 4-dimensional gauge theories. Utilizing 
the method, we show that pure-glue QCD exhibits very small positive dynamical tendency 
for self-duality. Speciflcally, using Wilson lattice deflnition of the theory and overlap-based 
deflnition of the fleld strength tensor, we flnd that the absolute X-distribution is mildly 



^In fact, as first shown in Ref. [15], the chiral polarization effects observed in works that followed Ref. [TT] 
and aimed at studying the merits of instanton liquid model of QCD vacuum [16] . are almost entirely due to 
kinematics. 
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convex with Ca ~ 0.018. This means that in head-to-head comparisons of samples chosen 
from QCD-generated field strength and from ensemble of its statistically independent dual- 
ity components, the probability for the former being more self-dual is only by Ca/2 ~ 0.009 
larger than equal chance. We also construct the effective low energy field strength based 
on overlap eigenmode expansion, which is found to have comparable degree of dynamical 
self-duality. 

2. The Absolute X— Distribution Method. We begin with a brief description of the ab- 
solute X-distribution method. Detailed discussion is given in Ref. [H]. Consider a quantity 
Q with values in some vector space that can be decomposed into a fixed pair of equivalent 
orthogonal subspaces, i.e. Q = Q1+Q2 with Qi ■ Q2 = 0. The dynamics of Q is governed by 
the probability distribution Vf{Qi, Q2) which is assumed to be symmetric with respect to the 
two subspaces, i.e. Vf{Qi,Q2) = 'Pf{Q2,Qi)- In field theory apphcations the distribution 
Vf{Qi, Q2) descends from the distribution of field configurations defined by the local action. 

We are interested in characterizing the property of polarization, namely the tendency 
for asymmetry in the contribution of the two subspaces in the preferred values of Q. The 
contribution of each subspace in sample Q is measured by magnitudes = \Qi\ of their 
components, and one can thus reduce the full distribution Vf{Qi, Q2) to the distribution of 
magnitudes Vb{qi,q2) for this purpose. Moreover, the asymmetry in the contribution of the 
two subspaces in Q can be assessed by the ratio q2/qi which can be cast into a normalized 
variable by considering a combination such as [H] 



X = - tan ^ f — ) - 1 = Xriqi, 92) x E [-1, 1] (1) 

TT \qiJ 

namely the reference polarization coordinate. Notice that x = — 1 if sample Q is strictly 
polarized in first direction, x = 1 if it is strictly polarized in second direction, and x = 
for strictly unpolarized sample qi = q2- The marginal distribution of reference polarization 
coordinate in distribution Vb{qi,q2), namely 

Pr{x) = j dqi j dq2Vb{qi,q2) 5(x - Xr{qi,q2)J (2) 

is the reference X- distribution [11] . It can be viewed as detailed but kinematical polarization 
characteristic of Vf{Qi,Q2) since its qualitative properties depend on the choice of X^. 

Rather than the distribution of reference polarization function (coordinate) Xr(x) = x, 
the absolute X -distribution Pa{.X) represents a distribution of the polarization function 

■Xa{x) = 2 r dyP:{y)-l (3) 
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where the reference X-distribution P"(x) is computed for the dynamics of statistically in- 
dependent components, namely 

/•oo 1*00 

K{Qi,Q2) = Piqi)piq2) P{q)= j dq2Vb{q,q2) = j dqiVb{qi,q) (4) 

Jo Jo 

There are three important points to note about Pa{X) so constructed, (i) Xa{x) is a 
valid polarization measure in exactly the same sense as x itself is: it quantifies the degree 
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of polarization in the sample, (ii) As one can easily check, the distribution of Xa{x) is 
uniform in uncorrclatcd distribution V^{qi, ^2) and thus the absolute X-distribution directly 
measures polarization tendencies relative to statistical independence. One can best see the 
differential nature of this direct comparison from the explicit form 

(Hi) The distribution Pa{X) is invariant under the choice of reference polarization function 
(coordinate) used in its construction, and is in this sense absolute. 

Taken together, properties (i-iii) indicate that the absolute X-distribution represents 
a dynamical polarization characteristics. This differential information can be reduced to 
integrated form by computing the associated correlation coefficient of polarization, namely 

Ca = 2 j dX \X\ PAiX) - 1 (6) 

The positive value of Ca G [—1, 1] indicates that the dynamics enhances polarization rela- 
tive to statistical independence, while the negative value implies its dynamical suppression. 
Indeed, the precise statistical meaning of Ca is as follows. Consider the experiment wherein 
samples are independently drawn from Vb{qi,q2) and V^{qi,q2) in order to be compared, 
head to head, with respect to their degree of polarization. The result of this experiment is 
the probability Ta that a sample drawn from ensemble in question is more polarized than 
sample drawn from ensemble of statistically independent components. Then Ta = (C^-|-l)/2 
thus implying the above interpretation of Ca- 

3. Duality Decomposition. Our object of interest is the Euclidean field-strength tensor 
F{x) of SU(N) gauge theory in four dimensions. To analyze its dynamical polarization 
properties with respect to duality, we decompose it into its self-dual and anti-self-dual 
parts, i.e. 

F^Fs + Fa with Fs = l{F + F) Fa = \{F-F) (7) 
where the dual tensor F is defined in a standard way, namely 

Fixv = \^iivpaFpa ei234 = 1 (8) 

The above split represents an orthogonal decomposition of the field-strength tensor relative 
to the scalar product 

A-B = \Y,^xA^,B^, (9) 

with A = {Afj_^}, B = {Bf^^} representing arbitrary antisymmetric collections oi N x N 
Hermitian matrices. Thus, to evaluate dynamical tendency for self-duality in SU(N) gauge 
theory we apply absolute polarization methods to objects {Qi,Q2) = {Ps,Pa), which in 
practice involves statistical analysis of pairs {qi, 52) ^ {\Fs\, \Fa\) with the norm defined by 
the scalar product above. 



PAiX) 
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4. Lattice QCD Setup. Our exploratory calculations will be performed in the context 
of pure glue SU(3) lattice gauge theory with Iwasaki action [17j. Table [1] summarizes the 
parameters of the ensembles. Lattice scale has been determined from the string tension 
following the methods and results of Ref. [H]. Ensembles E2, E3, Eg, E^ and Ej have the 
same physical volume, and serve to investigate the continuum limit. Ensemble E^ has larger 
physical volume to check for possible finite size effects. 

We will utilize the local lattice definition of gauge field strength tensor based on the 
relation [HI |20l [21] 

tr,a^,Do,o(f/(a)) = F^,{Q) + 0{a^) (10) 

with = c^{p,r) being generically non-zero for overlap Dirac operator D = D'^''^ [22]. The 
above expression is valid for discretization U{a) of arbitrary classical field A^, and tr^ denotes 
trace over spin indices only. Parameters r = 1 and p = 26/19 were used for calculations 
throughout this paper. Local values of Ffj,^{x) in equilibrium lattice QCD configurations 
were obtained by evaluating the needed matrix elements of D using point sources. 

One advantage of the above definition for the field strength is that it can be naturally 
eigenmode-expanded to provide an effective low-energy representation. To do that, we follow 
the same strategy as one used in definition of the effective topological density in Ref. [8] . In 
particular, rather than truncating the expression for lattice field implied by (JTOl) . we take 
advantage of the fact that tr cr^j^ = and rather expand 



-P^fiui^-^) — tr^ (7^1, Dx^x 



— tr, (2p - Dx.x) 
Ct 



Indeed, while the expansion of the former gives lowest weight to low-lying modes and highest 
weight to the modes at the cutoff, the latter reverses this and gives highest modes zero weight 
as naturally desired. The effective field at the fermionic scale A is then 



a,b 



i E '^^ 

|A|<Aa 



A) 



a 



(12) 



where ipx{x) is the eigenmode of D with eigenvalue A and a,b are color indices. 

Low-lying eigenmodes needed for the construction of effective field strength were com- 
puted using our own implementation of the implicitly restarted Arnoldi algorithm with 
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80 




39 


96 


19 




20 


20 



E2 
E3 
Eg 

^4 
Eq 



12^ 
16^ 
20^ 
244 
32^ 
32^ 



;i.32fm)4 
^1.32 fm)4 
;i.32fm)4 
;i.32fm)4 
;i.32fm)4 
;i.76fm)4 



0.110 fm 
0.0825 fm 
0.066 fm 
0.055 fm 
0.041 fm 
0.055 fm 



2.530 
2.725 
2.892 
3.0375 
3.278 
3.0375 



Table 1: Ensembles used in the overlap calculations of Ff^^{x). Column A^^^^gg shows the 
number of configurations for which local field strength (11 01) has been computed. Column 
^config lis^s number of configurations with computed eigenmodes when available. 
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Figure 1: Absolute X-distribution in duality and reference X-distribution of statistically 
independent duality components for ensemble Eq. 



deflation [23]. The zero modes and about 50 pairs of lowest near-zero modes with complex- 
conjugated eigenvalues were obtained for each configuration. In Table [1] we indicate the 
ensembles with eigenmodes available for these purposes. One point to note here is that, 
as one can easily check, the terms associated with chiral zero modes are strictly self-dual. 
While the relative weight of this contribution is zero both in the infinite volume limit at 
fixed lattice spacing, and in the continuum limit at fixed physical volume, the zeromodes do 
contribute in finite lattice calculations and are included. 

5. Main Results. We have computed absolute X-distributions in duality for all ensem- 
bles listed in Table [H The plot of Fig. [T] shows the result for ensemble Eq which has the 
largest physical volume. Note that the associated kinematical background represented by 
the reference X-distribution of statistically independent components is also shown0 As can 
be clearly seen, despite the completely "unpolarized-looking" kinematics, there is a slight 
dynamical tendency of the gauge field strength tensor to polarize itself. Indeed, the abso- 

^As emphasized in Ref. jM], the combination of absolute X-distribution and reference X-distribution 
of uncorrelated components provides all that is necessary to construct X-distribution in arbitrary reference 
frame (polarization function) . In that sense, it stores a complete information on polarization, dynamical and 
kinematical. 
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lute X-distribution shows a small excess of probability near the extremal values of absolute 
polarization coordinate. This observation represents the central message of this work since, 
as we discuss below, the data suggests that such dynamical behavior persists both in the 
continuum and infinite volume limits. 

Focusing on absolute X-distributions from now on, we show them more closely for all 
ensembles with fixed physical volume in Fig. |2l As can be directly inspected, these distri- 
butions are all mildly convex, implying small positive dynamical tendency for polarization 
in duality. Despite the fact that this tendency slowly decreases toward the continuum limit, 
the data strongly suggests that a finite positive continuum limit does exist. In Fig. [2] we 
also show a fit to quadratic polynomial in lattice spacing, which we take to facilitate this 
continuum extrapolation. 

Although the physical volume in the above continuum extrapolation is rather small, the 
finite volume corrections for observables at hand turn out to be negligible. Indeed, in Fig. [3] 
we show the comparison of absolute X-distributions for ensembles £'4 and Eq that have 
identical couplings (a = 0.055 fm) but with Eq representing a significantly larger physical 
volume. As can be seen in this close-up view, the two distributions are rather difficult to 
distinguish one from another. 

The above results make strong case for the following proposition that we suggest for 
further investigation. 

Proposition 1: Absolute X -distributions for duality in SU(3) pure glue lattice gauge theory 
are convex. This property, and thus the associated positive tendency for (anti) self -duality, 
will persist in the continuum limit. 

It needs to be emphasized that the observed polarization tendency, as objectively quanti- 
fied by the associated correlation coefficient, is very small. In particular, Ca ~ 0.018 for 
continuum limit taken using Iwasaki lattice gauge action and overlap-based definition of the 
field-strength tensor. 

6. Effective Field Strength Tensor. It is revealing to compare dynamical polarization 
tendencies of fully fluctuating field strength to those of effective fields defined by equation 
f|T2|) . In FigJH (top right) we show the comparison of absolute X-distributions for F and 
F^, A = 1000 MeV, in ensemble i?4. The construction of effective field strength in this case 
involved the inclusion of 47 overlap near-zero modes on average. As one can see, while the 
dynamical polarization tendency has somewhat increased in the effective field, the two cases 
do not differ qualitatively at all. This happens despite the fact that the associated reference 
X-distributions, which are only kinematic, differ significantly (top left of FigJl]). In the 
bottom panel of FigJUwe added to the continuum extrapolation plot of Figj2]the correlation 
coefficients Ca at A = 1000 MeV for ensembles with available eigenmodes. This shows that 
the correlation coefficients in the effective field are of comparable magnitudes to those for 
full field strength, and more so in the continuum limit. 

Effective fields at fermionic scale A suppress short distance fluctuations regardless of 
whether they contribute dynamically or represent pure noise. The above comparison thus 
speaks to the fact that absolute polarization measure has little sensitivity to the presence of 
noise, as emphasized in the opening remarks, and to the fact that the dynamical self-duality 
effect is mainly due to low-energy fields. 
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Figure 2: Absolute duality X-distributions for all ensembles at fixed physical volume V = 
(1.32 fm)*^. Lattice spacing is decreasing in lexicographic order of plots. The last panel 
shows the corresponding lattice spacing dependence of the polarization correlation coefficient, 
together with the quadratic fit for extrapolation to the continuum limit. All data points in 
the above plots have errorbars too small to be clearly seen. 
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Figure 3: The comparison of absolute X-distributions at the same lattice spacing but dif- 
ferent volumes {E^ and £"6). 



7. Discussion. In this work we constructed dynamical measures of self-duality in the gauge 
field. These measures are based on the notion of absolute X-distribution [1^ which in turn 
represents a differential correlational characteristic of polarization. From the conceptual 
standpoint, this is an important step forward from initial approaches of Refs. pTl [13] that 
are only kinematical, and thus arbitrary. Indeed, not even a sign of dynamical tendency can 
be deduced from such results. 

Using the above framework, we performed a detailed study of self-duality in pure-glue 
SU(3) gauge theory using lattice regularization. Our results indicate that this theory in- 
volves positive dynamical tendency for self-duality, meaning that QCD dynamics produces 
enhancement of self-duality relative to statistical independence in duality components. How- 
ever, the observed effect is very weak. Indeed, the measured correlation coefficient is only 
Ca ~ 0.018 while, for example, it is straightforward to prescribe dynamics for measured 
marginal distributions of duality components, that would enhance this correlation at least 
20-30 times. It is significant in this regard that the same is true for low-energy effective 
fields obtained via eigenmode expansion of the field-strength tensor, despite of the associated 
"noise reduction" . We are thus led to conclude that self-duality, and hence classicality, does 
not manifest itself as a significant feature of pure-glue QCD dynamics. This is consistent 
with previous arguments as well as with lattice QCD studies [HI El [251 ID] indicating 

an intrinsically non-classical paradigm of QCD topological charge fluctuations. 
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Figure 4: Comparison of dynamical polarization properties for field strength tensor F and 
the effective field strength F^, A = 1000 MeV. See discussion in the text. 
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